On the scaling laws for the largest Lyapunov exponent 
in long-range systems: A random matrix approach 



O 

o 

(N 
D 



Celia Anteneodo* and Raul O. Vallejos^ 
Centra Brasileiro de Pesquisas Fisicas, R. Dr. Xavier Sigaud 150, 
22290-180, Rio de Janeiro, Brazil 

We investigate the laws that rule the behavior of the largest Lyapunov exponent (LLE) in many 
particle systems with long range interactions. We consider as a representative system the so-called 
Hamiltonian q-XY model where the adjustable parameter a controls the range of the interactions of 
A'' ferromagnetic spins in a lattice of dimension d. In previous work the dependence of the LLE with 
the system size A^, for sufficiently high energies, was established through numerical simulations. In 
the thermodynamic limit, the LLE becomes constant for a> d whereas it decays as an inverse power 
law of A'^ for q < d. A recent theoretical calculation based on a geometrization of the dynamics 
is consistent with these numerical results. Here we show that the scaling behavior can also be 
explained by a random matrix approach, in which the tangent mappings that define the Lyapunov 
exponents are modeled by random simplectic matrices drawn from a suitable ensemble. 
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I. INTRODUCTION 

Dynamical systems of many particles interacting via 
long range forces can exhibit interesting anomalies such 
as super-diffusion ||l|-§l, metastable states Q and non 
Boltzmann-Gibbs distribution functions Q. Special in- 
terest in such systems has arisen recently in connection 
with the non-extensive generalization of statistical me- 
chanics introduced by Tsallis 

The existence of a dynamics makes the systems men- 
tioned above very attractive because it is possible, in 
principle, to associate the properties of the thermody- 
namic states with features of the many particle phase 
space. As a remarkable example, let us mention the 
"Topological Hypothesis" , which relates thermodynamic 
phase transitions to topological changes in the structure 
of phase space . 

A dynamical model with an adjustable interaction 
range a, and allowing extensive numerical and analyti- 
cal exploration has been recently introduced jsj-p^ . The 
model consists in a periodical c?-dimensional lattice of TV 
interacting rotators moving on parallel planes. Each ro- 
tator is restricted to the unit circle and therefore it is 
fully described by an angle < 6*^ < 27r and its conjugate 
momentum L^, with i = 1, . . . , N. The dynamics of the 
system is governed by the Hamiltonian 
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where the coupling constant is J > and, without loss 
of generality, unitary moments of inertia are chosen for 



all the particles. Here Vij measures the minimal distance 
between the rotators located at the lattice sites i and j. 
The Hamiltonian (|^) describes a classical inertial XY fer- 
romagnet. It contains as particular cases the mean- field 
version (a/d — 0) and the first-neighbors case, recovered 
in the a/d — > oo limit. 

The systems a-XY characterized by interaction ranges 
< a/d < 1 do not have a well defined thermodynamic 
limit, e.g., the specific (per particle) energy diverges when 
N ~> oo. A proper thermodynamic limit is defined by in- 
troducing a scaling parameter N Jill, which depends on 
N, a, and d MHI 
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In the large N limit one has 

jyl-a/d 0<a <d 

N{a/d)r^{ In TV a^d 

Q{a/d) a > d 



(3) 



with O a function of the ratio a/d only. Specific energy- 
like quantities must be rescaled by N. At the dynamic 
level, time, hence inverse Lyapunov exponents, have to 
be scaled by TV^/^ 

A completely equivalent description is obtained by 
working with the already scaled Hamiltonian 
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This kind of scaling of the strength of the interactions, 
common in standard mean-field discussions, has been ap- 
plied to the study of the mean-field case of the present 
model in Refs. Since the Hamiltonian leads to 

the same results as (|l|), but avoiding further rescalings, 
all our considerations from here on will be related to the 
already scaled Hamiltonian (^. 

The a-XY ferromagnet has been subject of several nu- 
merical and analytical studies. References and 
P,p^ are dedicated to the cases d = 1 and d — 2,3, re- 
spectively. The mean-field problem is discussed in Refs. 
p^ , p^ while the opposite limit of first-neighbor inter- 
actions and d — 1 can be found in For long- 
ranged interactions, a/d < 1, the system displays a sec- 
ond order phase transition from a ferromagnetic state to 
a paramagnetic one at a certain critical specific energy 
£c (ec — 0.75 J) |l5| , [l6|| . It is hkely that systems with 
1 < a/d < 2 also undergo a second order phase transi- 
tion but the critical energy may depend on a/d, as has 
been shown numerically for d = 1 [ p^ . For a/d > 2, 
all systems behave similarly to the first-neighbor model, 
where there is no order nor phase transition for finite 
energies in the thermodynamic limit PJT^]. 

Here we are concerned with the high energy phase, i.e., 
with energies above Ec- In this disordered regime, the ki- 
netic energy is much larger than the bounded potential 
energy, the rotators are weakly coupled and the dynamics 
is weakly chaotic. In Refs. [p|-p^ largest Lyapunov expo- 
nents (LLE) were calculated numerically. It was found 
that, in the thermodynamic limit, and for large enough 
energies, the largest Lyapunov exponent remains positive 
and finite for short-range interactions {a/d > 1) but van- 
ishes as an inverse power law of the system size in the 
long-range case < a/c? < 1. 

Recently Pettini and co-workers developed a theoret- 
ical method which, in principle, allows to obtain the 
scaling behavior of the LLE analytically (see for a 
review). In this approach the phase space trajectories 
are mapped onto a geodesic flow in configuration space 
(equipped with a suitable metric). It is assumed that the 
curvature along an ergodic geodesic can be modeled as a 
Gaussian process. Then the LLE is expressed in terms 
of the mean and variance of the process. These parame- 
ters are calculated as microcanonical averages of suitable 
dynamical functions. There are several works where the 
method was applied to the a-XY model. The scaling be- 
havior of the LLE in the extreme cases a — > oo and a = 
was found in Refs. and respectively. Very re- 
cently, Firpo and Ruffo succeeded in calculating the 
LLE scaling laws for any interaction range Q < a/ d < \. 

It is our purpose here to present a simple alternative 
procedure, based on a random matrix formulation, which 
allows to derive the dependence of the LLE on the size 
of the system N , the range of the interactions a and 
the lattice dimension d. This procedure is based on the 
ideas introduced by Benettin pTI] in the discussion of two 



dimensional billiards, and later extended to interacting 
many particle systems p2|,p3[. 



II. THE LARGEST LYAPUNOV EXPONENT 

The equations of motion generated by the Hamiltonian 
0) are 
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for i = 1, . . . ,N. Discretizing the time axis into steps 
At one obtains the stroboscopic map relating angles and 
momenta at successive discrete times 



9' = + L, At 
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For the purpose of discussing Lyapunov exponents one 
has to consider the tangent map T, i.e., the linearized 
version of Eqs. (0), given by 



9' = 59, + 5U At 
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In matrix form these equations read 
f 5e'\ ft lAt \ f 50 
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where all submatrices are of size iV x iV, 1 being the 
identity. The matrix b is given by 



cos{e[-er)/Tf^ for 3^1 



N 



J2 cos(0:-0^)/,f, for j^i 



and the perturbation parameter e is 
e ^ J/ N{a/d) . 
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For long-range interactions it is clear that e goes to zero 
in the thermodynamic limit. In the short-range case one 
has e ~ J. In order to treat e as a small parameter in 
both cases, we will consider the limit J ^ when neces- 
sary. 

The LLE can be defined by the limiting procedure 



Amaa: = 1™ ^-ln||T^|| , 



(11) 
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with ^ an arbitrary vector and T = T„T„_i . . . T2T1 is 
the product of n tangent maps calculated at successive 
points of the discretized trajectory. Using the euclidean 
norm, Eq. (O) is rewritten as 
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n^oo 2nA.t 
the superscript t indicating "transposed" . 

III. RANDOM MATRIX APPROACH 

The random matrix approach is based in the model- 
ing of the tangent mappings Tk by a sequence of non- 
correlated random simplectic matrices mimicking the es- 
sential properties of the chaotic dynamics. In the stan- 
dard procedure |^ one replaces the short-time Tk by 
finite time random matrices Rk having the same struc- 
ture 
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In previous treatments the time scale r has always been 
ignored by setting it to one. Here we prefer to keep track 
of r, as later we will argue that it is related to energy. 
It is fixed as follows. The time scale r must be chosen 
small enough in order to preserve the short time structure 
of the tangent maps. However, it cannot be too small, 
as we assume that consecutive tangent maps are statis- 
tically independent. Thus t must be an intermediate 
scale, of the order of the correlation time for the tangent 
maps. (This time scale is analogous to the correlation 
time of the Gaussian process modeling the fluctuations 
of the curvature in the geometric method [l^.) 

The symmetric matrices are the random analogs of 
hk and are assumed to be independent, i.e., (a„am) = 
{an)(am), unless n = m. Except for the symmetry re- 
striction, the elements of a given matrix ak are inde- 
pendent |p5| . 

The probability distributions of the elements are 
obtained from bij (||) by considering that the angles ■ are 
independent from each other and uniformly distributed 
in [0,27r). This assumption, reasonable in the high en- 
ergy phase, implies that: (i) The average of each is 
null, (ii) The information about the range of the inter- 
actions is embodied in the variance of each a^j, which 
depends on the distance between sites: 

1 



{iau)') = lN{2a/d) 
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In this way one has defined a crude although non triv- 
ial statistical model whose validity has been shown in 
previous works pl|-p3[. 



As a consequence of the assumptions made one has the 
property (afefl^) = 7I, which will be useful for evaluating 
averages. In our case 
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Within this model the expression for the LLE is obtained 
by averaging over different realizations of sequences of 
random matrices 



Xmax = lim -^{\n£_*Rl 
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Assuming that the distribution of LLE's over the ensem- 
ble of sequences is narrow, one can interchange the av- 
erage and the logarithm. Then the averaging scheme is 
reduced to a sequence of averages over each matrix dis- 
tribution 



lim T^nei- 
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These averages have already been calculated by Parisi 



and Vulpiani |23| (see also |24|). Instead of just recalling 
their results, we prefer to exhibit a different derivation, 
which not only makes the paper self-contained, but may 
be interesting by itself. Notice that the first average can 
be immediately done, the result being the symmetric ma- 
trix 
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with 
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The remaining n — 1 averages can be done iteratively, ob- 
taining at each step a symmetric matrix with the same 
structure as (llq). The final result is 



{R{{...)R,) = 



ant 
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where the coefficients v, fJ,,(T are calculated by means of 
the recurrence relation 



(21) 



Then it can be easily checked that the LLE is related to 
the maximum eigenvalue Lmax of the 3x3 matrix above 
through the formula 
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After solving the cubic eigenvalue equation we expand 
0, 



so that 



= 1 + (27eV)i/3 + ... , 



-(27eV)i/3 



Finally, substituting e and 7 by their definitions (|1_ 
(plsj), respectively, one gets the compact expression 



A oc 72/3^1/3 
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(24) 
and 

(25) 



However, for the purpose of comparison with previous 
works, it is convenient to explicit the dependence. Re- 
calling the asymptotic expression (|^) for TV, we arrive 
at 



l/Ari/3 
(lniV/7V)i/3 

l/(lniV)2/3 
constant 



0<a/d< 1/2 

a/d^ 1/2 
1/2 <a/d<l . 
a/d = 1 
1 < a/d 

(26) 



This scaling law for the LLE can also be written as 
Knax ^ l/A^" with 



1/3 
2{1 ~a/d)/3 




0<a/d< 1/2 
1/2 <a/d<l 
1 < a/d 
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the case a/d = I being marginal. 

We expect the scaling above not to depend on the de- 
tails of the dynamics, i.e., it should be typical of systems 
with couplings of the form 1/r", e.g. classical n- vector 
ferromagnets (of which n = 2 is the present case), as 
long as the perturbation d has zero mean. Systems for 
which the average perturbation is non-zero belong to a 
different universality class, and alternative scaling laws 
are expected pi|- 



IV. CONCLUDING REMARKS 

The scaling behavior of Xmax with N 00 and J — > 
[Eq. (p6|)] is exactly the same as that obtained by using 
the geometric method The agreement can be also 

extended to the energy domain (e — > 00) by relating the 
time scale r to energy. Given that the potential energy is 
bounded, when e — > 00, the total energy and the kinetic 
energy are essentially the same. In this regime, chang- 
ing the time scale is equivalent to a change in the kinetic 



energy, so that we have r oc e"^/^. Thus we obtain the 
scaling law Xmax oc e~^l^ . 

The theoretical results ( |27| ) agree with the numerical 
ones obtained numerically in Refs. [p|-p^. There are some 
deviations which are consistent with finite size effects, as 
argued in However, one should not discard the pos- 
sibility that the scaling laws are not exactly those derived 
in this paper (or in |20|). The differences with numeri- 



cal calculations might be due to the fact that both the 
geometric and the random matrix approaches assume er- 
godicity and fast (exponential) decay of the correlations. 
We do not know at present if the dynamical system fully 
satisfies these hypotheses. Eventually this issue will be 
decided when simulations on larger systems are available. 
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